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^sD ! Abstract 

We consider brane world models, which can be constructed in the five-dimensional 
fNj ■ Brans-Dicke theory with bulk scalar field potentials suggested by the supergravity theory. 

, For different choices of the potentials and parameters we get: (i) an unstabilized model 

with the Randall-Sundrum solution for the metric and constant solution for the scalar 

■ field; (ii) models with flat background and tension-full branes; (iii) stabilized brane world 
Q\ . models, one of which reproduces the Randall-Sundrum solution for the metric and gives 

an exponential solution for the scalar field. We also discuss the relationship between 

■ solutions in different frames - with non-minimal and minimal coupling of the scalar field. 

> 

CO , 

1 Introduction 

CN 

Brane world models and their phenomenology have been widely discussed in the last years 
[U [21 [3J HI [5J [6]. One of the most interesting brane world models is the Randall-Sundrum 
model with two branes, - the RSI model [TJ. This model solves the hierarchy problem due 
to the warp factor in the metric and predicts an interesting new physics in the TeV range of 
Ph! energies. A flaw of the RSI model is the presence of a massless scalar mode, - the radion, which 
describes fluctuations of the branes with respect to each other. As a consequence, one gets a 
scalar-tensor theory of gravity on the branes, the scalar component being described by the 
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radion. The radion coupling to matter on the negative tension brane contradicts the existing 
restrictions on the scalar component of the gravitational interaction (see, for example, review 
[8]), and in order the model be phenomenologically acceptable the radion must acquire a mass. 
The latter is equivalent to the stabilization of the brane separation distance, i.e. it must be 
defined by the model parameters. The models, where the interbrane distance is fixed in this 
way, are called stabilized models, unlike the unstabilized models, where the interbrane distance 
can be arbitrary. 

Such a stabilization of the interbrane distance can be achieved, for example, by introducing 
a five-dimensional scalar field with brane potentials [9] . A disadvantage of this approach is that 
the backreaction of the scalar field on the background metric is not taken into account. This 
problem is solved in the well-known model proposed in [10]. However the background metric of 
this solution is rather complicated and differs significantly from the Randall-Sundrum solution. 
An interesting problem is, whether it is possible to find a stabilized model, where the Randall- 
Sundrum form of the metric is retained despite the interaction with the scalar field. 

A solution to this problem based on a non-minimal coupling of the scalar field to gravity 
was put forward in [llj . where the form of the metric was found to be the Randall-Sundrum 
one, whereas the solution for the scalar field again turned out to be rather complicated. 
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One of the standard forms of the non-minimal coupling of the scalar field to gravity is 
the linear coupling to scalar curvature used in the Brans-Dicke theory of gravity. Brans-Dicke 
theory in the brane world context was already discussed in the literature, see, for example, 

nana mi. 

In the present paper we use this type of coupling to construct a stabilized model with the 
Randall- Sundrum solution for the metric, the solution for the scalar field being also a simple 
exponential function. Our paper is organized as follows. First, we present a method for con- 
structing different background solutions in five-dimensional Brans-Dicke theory by considering 
bulk and brane scalar field potentials of a special form and examine their correspondence with 
the solutions in the Einstein frame. In particular, we show that the solutions with these special 
scalar field potentials in the Jordan frame correspond to the solutions in the Einstein frame 
with potentials suggested by the supergravity theory. Then we discuss models which can be 
obtained with different choices of the potentials and parameters. And finally, we discuss the 
obtained results. 



2 Setup 

Let us denote the coordinates in five- dimensional space-time E = M4 x S l /Z2 by {x M } = 
y}, M = 0, 1, 2, 3, 4, \i = 0, 1, 2, 3, the coordinate x A = y, —L < y < L parameterizing the 
fifth dimension. It forms the orbifold, which is realized as the circle of the circumference 2L 
with the points y and — y identified. Correspondingly, the metric qmn and the scalar field 
satisfy the orbifold symmetry conditions 

g^(x, -y) = g^(x, y), g^x, -y) = -g^(x, y), (1) 
gu(x, -y) = gu(x, y), <P(x, -y) = <p(x, y). 

The branes are located at the fixed points of the orbifold, y = and y = L. 
The action of the brane world models can be written as 



S = d x dyy^g 



L 



^■g^du^-V^) 



(2) 



g\i((f))d x — y — g\2(4>)d x. 
y=o Jy=L 



Here V(<p) is a bulk scalar field potential and Xi y 2{<P) are brane scalar field potentials, g = 
detg^v, g^ v denotes the metric induced on the branes and u is the five-dimensional Brans-Dicke 
parameter. The formal restriction on this parameter in the five-dimensional Brans-Dicke theory 
is uj > — I (which can be easily seen, for example, from the formulas of paper [15]), and further 
restrictions on its value can be obtained only after studying an effective four-dimensional theory 
for a certain background solution. The only difference from the classical Brans-Dicke theory is 
the presence of the bulk scalar field potential V(<f)) and the branes. The signature of the metric 
g M N is chosen to be (-, +, +, +, +)• 

The standard ansatz for the metric and the scalar field, which preserves the Poincare in- 
variance in any four- dimensional subspace y = const, looks like 

ds 2 = e^rj^dx^dx" + dy 2 , (f>(x, y) = <j>(y), (3) 
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T]^ denoting the flat Minkowski metric. If one substitutes this ansatz into the equations 
corresponding to action one gets a rather complicated system of nonlinear differential 
equations for functions a(y),4>(y): 



3a" 6 + 
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L'\2 



h" - a'<f>' + ^8(y) + ^S(y - L) = 0, 
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V) + 4a 



0. 



)" + 4 ( xV)-4 ( T"-lO(a') 2 -^(0 / ) 2 



(4) 

(5) 
(6) 



IdV 1 dXi . 1 d\ 2 _ N 

where prime denotes the derivative with respect to extra dimension coordinate y. 
We will consider a special class of bulk potentials, namely 



V{<f>) = -(3cu + 4) 



(4cu + 5) 0F 2 + 20 2 F 



dF 



(dF 



(7) 



where F = is a function. This structure of the potentials is similar to the one introduced 
in [16] for the case of the minimal coupling of the scalar field to gravity. Although the potential 
in its general form ([7j) looks rather complicated, the resulting potentials for different choices of 
the function F, as we will see below, look quite natural. 

One can check that in this case any solution of the equations 



0F - 3(f) 2 



,dF 



a 



(oj + 1)F + 



dF 



(9) 



also satisfies (J4]), ((Sj) and (jSJ) in the interval [0, L], provided that the following boundary condi- 
tions on the branes are satisfied 



(3oj + 4) F| 



Ai 



y=+o 



(10) 



(3w + 4)^1^-0 = 71 



x dF. 
(3^ + 4) — — |y =+0 



l rf(W<ft) 
4 d6 



dF l djXM 



'ID 



(12) 



(13) 



It is necessary to note that the symmetry conditions ([I]) were used to obtain these relations. 

Thus we get first order differential equations instead of the initial second order differential 
equations. The situation is analogous to that in papers pUl EE] , where the bulk and the brane 



3 



potentials for the scalar field minimally coupled to five- dimensional gravity are chosen in an 
appropriate way. 

It is a common knowledge that action ([2]), called the action in the Jordan frame, can be 
brought by a conformal transformation to an action with the scalar field minimally coupled to 
gravity, which is called the action in the Einstein frame. Indeed, if in action (j2J) we put <ft = 
2M 3 exp(— ^) with a = 6M 3 / 2 (& + 4:/3) 1 / 2 and make a conformal rescaling g^N ~^ eyL v(^)9MN 
allowing one to pass from the Jordan frame to the Einstein frame, we get a model, describing 
scalar field p minimally coupled to five-dimensional gravity, with the action 

S = j d 4 x ^ dy f 2M 3 i2 - l -g MN d M pd NP - V{p)^j - (14) 

- / \/-^Ai (p)d A x - / ^§X 2 (p)d i X , 

Jy=0 Jy=L 

where 

V(p) = e^(2M 3 e"^), \i(j>) = e^Ai(2M 3 e-^), i = 1, 2. 
Thus, if we a have a solution in the Jordan frame 

ds 2 = e 2 ^ Vtlu dx^dx u + dy 2 , <f)(x, y) = <f>(y), 

generated by potential (JTj), it can be transformed to a solution of the corresponding theory (fl4|) 
in the Einstein frame 

ds 2 = e^)-^\ vdx ^ + e^dy 2 , p(y) = In (M\ . 



To bring the metric of this solution to the standard form ([3]), we have to pass to a new coordinate 
z of th 
frame 



z of the extra dimension according to dz = exp(—^-)dy. Thus, we get a solution in the Einstein 



ds 2 = e^M-^rj^dx^dx" + dz 2 , p(z) = p{y{z)), 
in theory ffl4|) with potentials 



dW x 2 



V(p(z)) = - ^r- —,W 2 , (15) 

^ v ;; 8\dp J 24M 3 ' v ; 

W(p(z)) = -^- e £ ^F(2M 3 e-'^), (16) 
3 

\{p{z)) = ±W(p(z))\ z=Zi , i = l,2, (17) 

which have exactly the same form as the potentials of the model discussed in [TUl [IS] , where 
their form was suggested by the supergravity theory. Thus, the method of finding solutions in 
the Brans-Dicke theory is equivalent to the one discussed in these papers, which is quite clear, 
because in both cases the second order differential equations are reduced to the first order ones. 
Nevertheless, if we study a Brans-Dicke theory, it is more convenient to look for solutions in 
the Jordan frame than to transform them from the Einstein frame. 

Another point is that solutions in one frame (Jordan or Einstein) may be elegant and 
interesting in view of the hierarchy problem, whereas in the other frame it is not the case. For 
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example, if in the Jordan frame we take polynomial functions F, in the Einstein frame we get 
solutions for theories with complicated potentials being a sum of several exponential functions. 
Of course, the converse is also true. For this reason it is also convenient to have a method of 
constructing brane world solutions directly in the Jordan frame. 

One more important point is that though solutions in different frames are equivalent, the 
corresponding theories in different frame are not equivalent, because if in the Brans-Dicke 
theory we had the minimal coupling of gravity to matter on the branes, then after the conformal 
transformation the scalar field would enter the term describing the interaction with matter on 
the branes (the so-called conformal ambiguity [17J). At the same time, since the model described 
in [10] is stable against fluctuations of gravitational and scalar fields, i.e. it is tachyon- and 
ghost-free (see [H]), the model described by ([2]) and © should possess this property too, 
although the bulk potentials for certain choices of F are unbounded from below. 

3 Specific examples 

3.1 The Randall- Sundrum solution 

Let us choose the function F{4>) in the following form 

F = B<fi, (18) 

where B is a constant. It follows from (jSj) that = const in this case. Bulk and brane potentials 
are 

V(0) = -^ 2 0l(3a; + 4) 2 , (19) 
Ai, 2 = +45(3cu + 4)03, (20) 

and solution for the warp factor is 

B i 

a> = — (3u + 4) 03 (21) 
3 

in the interval [0, L]. Making redefinition 

^B (3u + 4) 03 = -k, (22) 

= 2M 3 , (23) 

we get 

V((f>) = -12k 2 M 3 , (24) 

Ai, 2 = ±12fcM 3 , (25) 

which formally coincide with the original RSI solution [7]. Nevertheless the linearized theory 
in this background differs from that in the RSI model, because the fluctuations of the scalar 
field add an extra degree of freedom. In this case the model cannot be stabilized (i.e. the size 
of the extra dimension cannot be fixed) by adding positively defined potentials on the branes, 
since the solution for the scalar field does not depend on y (this point is discussed in detail in 
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subsection 3.3). Thus, this example is not interesting from the physical point of view, but it 
shows how the general method for constructing solutions works. 

It is not difficult to see that action fl2]) with potentials of the form (TT91) and (}2T)j) can be 
brought by a conformal transformation of the metric and coordinate transformations, described 
in Section 2, to the action of the unstabilized RSI model with minimally coupled scalar field. 
The only significant difference is that if in the Brans-Dicke theory we had the minimal coupling 
of gravity to matter on the branes, then after the transformation the scalar field would enter 
the term describing the interaction with matter on the branes. 



3.2 "Consistent" ADD scenario 

Let us consider the case, where u = 0, i.e. the kinetic term for the Brans-Dicke scalar field is 
absent. Let us also suppose that F = B/<f>. It is not difficult to check that 

V(<f>) = 0, (26) 

a' = (27) 

in this case. Equation (12?]) means that we have the flat five-dimensional background metric. 
Solution for 0, following from (jSj), has the form 

4> = ABy + D (28) 

in the interval [0, L], where D is a constant. From the boundary conditions one gets 

Ai j2 = TA, (29) 

where A is a constant defining the brane tensions. Finally we get 

, A, 
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\y\ + D. (30) 



Thus, we get the model with the flat five-dimensional background metric and tension-full branes, 
which was discussed in detail in paper [19]. This model can be easily stabilized by the same 
method as the one, which will be discussed in the next subsection. 

3.3 Stabilized brane world with the Randall- Sundrum solution for 
the metric 

Now let us consider the case F = const, i.e. it does not depend on the field 0. Bulk and brane 
potentials, corresponding to such a choice of F(4>), can be chosen to be (see ((Tj) and (|T0|) - (|T3|) ) 

V(cf>) = A0, (31) 

A lj2 = ±\<j>. (32) 
Let us suppose that A > 0. It is not difficult to check that from (j8"j)- (JT3"j) follows 

a = -k\y\, (33) 
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= Ce- U \ y \ (34) 

with 

U = \f (3co + AH4u + 5j' (35) 
fc=(w + l)u, (36) 

and the fine-tuning relation 




A = 4V-AW- — — . (37) 
V 4cu + 5 v 7 

Constant C is not defined by the equations. One can see that in the limit w^oowe arrive at 
the standard Randall-Sundrum solution. 

Now let us discuss stabilization mechanism which can be utilized in the case under consid- 
eration. We will follow the way proposed in [TO] and add stabilizing quadratic potentials on 
the branes, namely 

AAi, 2 = 71,2 (0 - vi :2 f , 71,2 > 0. (38) 
Such an addition will not affect equations of motion provided 

(fi\ y =o = v x , <p\ y =L = v 2 . (39) 

Thus, now the constant C appears to be defined and is equal C = v\, whereas the size of extra 
dimension is now defined by the relation 

u \V 2 J 

which is the same as in the model of paper [TO]. Obviously, this stabilization mechanism works 
only for solutions with the scalar field depending on the extra coordinate, in particular, for the 
solution in subsection 3.2, and does not work for the solution in subsection 3.1. We would like 
to emphasize that this mechanism differs somewhat from the Goldberger-Wise mechanism [9], 
because it takes into account the backreaction of the scalar field on the metric, and the size of 
the extra dimension is fixed by the boundary values of the former. 

Now let us find the relationship between the four-dimensional Planck mass and the param- 
eters of the theory. We assume that the brane at y = L is "our" brane. To this end one should 
choose a = —k\y\ + kL to make four-dimensional coordinates {x^} Galilean on this brane (this 
problem was discussed in detail in [20]). Naive considerations suggest (see, for example, [19J) 
that the wave function of the massless four-dimensional tensor graviton has the same form as 
that in the unstabilized Randall-Sundrum model, namely h® (x, y) = e 2 ° a^ u {x). Moreover, 
the form of the residual gauge transformation, which are left after imposing the gauge on the 
fluctuations of metric (see [T8~| [20]). also suggests the same form of the wave function. Thus, 
substituting the following ansatz 

9^u(x, y) = e 2a g^(x) } </>(x, y) = <j)(y), (41) 
g M (x,y) = l, g^(x,y) = 
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into action (J2J), we get 

S = J J>e tc dy J R {4) ^^d 4 x (42) 

and 

2M|, = f L <j>e*°dy = (e* kL - e~ uL ) . (43) 

J-L IK + U 

If uL < 1 we get a formula analogous to that in the unstabilized Randall-Sundrum model (see 
J6j[20]). To solve the hierarchy problem one needs kL to be of the order kL ~ 30. If one chooses 
relatively large u (for example, uj > 30), then k would go to the value corresponding to the 
unstabilized Randall-Sundrum model, namely 



*«v-£ (44) 

(compare with ( 1231) and ( 1241 ) whereas uL could be made less than unity (uL < 1), since 
u = k/{u + 1). Under this assumption the parameter v\ can be chosen to be of the same 
order as t>2- In this case the parameters of the model, made dimensionless with the help of a 
fundamental scale in the TeV range, do not contain a hierarchical difference. The situation 
turns out to be completely analogous to that in the model proposed in PTJ] . At the same time, 
solution for the warp factor in stabilized brane world model found in [TU] is quite complicated, 
which impedes the analysis of the equations of motion for linearized gravity (approximate 
solutions can be found only under certain assumptions and simplifications, see [18]). As for our 
case, one can think that though the general structure of action ([2]) is more complicated than 
that of the action used in pi)] , the simplicity of solutions ( 1331) and (1341) could result in simpler 
equations of motion for linearized gravity. 

Another advantage of the solution presented above is that one can use all the results, 
obtained for the case of the universal extra dimensions in the Randall-Sundrum model (i.e. 
if one allows additional fields to propagate in the bulk, see, for example, [21J and references 
therein), in our case too. This happens because of the equivalence of solutions for the warp 
factors in both models. Of course, it is true in the case of the standard coupling of five- 
dimensional gravity to matter in the bulk. But since the size of extra dimension in our model 
appears to be stabilized, one can think that this would allow us to avoid possible problems 
caused by the radion field, which are inherent in the unstabilized Randall-Sundrum model. 

Quite an interesting situation arises, if one chooses uj = —1. Although the Brans-Dicke 
parameter is negative, the model is stable, since uj > — |. Formulas (|331) - (|3T|) with u = — 1 
take the form 

a = 0, (k — 0), (45) 
which means that the five-dimensional background metric is flat, 

(f) = Ce~ uM , (46) 

u = y/^A, (47) 

A = 4 V /= A. (48) 

Thus we get a model, which is similar, to some extent, to that discussed in subsection 3.2. 
In order to have the hierarchy problem solved (in the way proposed in [19]) in case of TeV 
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range of fundamental five-dimensional physics, one should choose uL ~ 30, as in the Randall- 
Sundrum model, and C ~ e uL . The flaw of the case u = — 1 is that there appears a new 
hierarchy between v\ and v-i. Nevertheless, this choice of parameters can be interesting from 
the pedagogical point of view since it demonstrates another scenario with flat background and 
tension-full branes (and nonempty bulk). 

The solutions of this section can be easily related to solutions in the Einstein frame. For 
these solutions F = const, and V(p(z)) (see ( |T5l) and ( |T6l) ) is the Liouville potential. Such 
" dilatonic" brane worlds were widely discussed in the literature, see, for example (22J E3] . As 
we have shown, our solutions can be brought to the form of the general solutions of paper [22J 
by a conformal rescaling of the metric and an appropriate transformation of the extra dimension 
coordinate, the latter being necessary for retaining the same ansatz for the metric (of the form 

©)• 

Our results also give a simple solution to the problem of finding a stabilized brane world 
model with the Randall-Sundrum form of the metric, which was discussed in [UJ. In this 
paper a similar solution in a theory with the scalar field non-minimally coupled to gravity, 
which preserves the Randall-Sundrum form of the metric, was found. But this solution and 
ours cannot be transformed to each other by a redefinition of fields, rescaling of metric and 
coordinate transformations. 

We would like to note once again that the theories obtained by a rescaling of the metric 
are not equivalent, if one considers the standard coupling of gravity to matter on the branes, 
as we have already mentioned in the end of Section 2. Since in stabilized brane world models 
the fluctuations of the scalar field describe also the radion (see [18J), this ambiguity modifies, 
in particular, the radion coupling to matter on the branes. Because we do not know, which 
frame is the "real" one, there are no strong objections against choosing the Jordan one. In this 
connection, an interesting problem is to compare the physical consequences of the models in 
different frames, which can be transformed one into another in the absence of matter on the 
branes. 



3.4 Power law solutions 

Let us consider the case 

F = B0 n , (49) 

where n ^ | and ti^O (these cases were discussed above). The corresponding bulk potential 
has the form 

V{4>) = -(3uo + A)B 2 <p 2n+l [4w + 5 + 2n - 3n 2 ] , (50) 

whereas the fine-tuned brane potentials can be easily obtained from (fTDl - ([T3"l) . Solutions for 
the scalar field and the warp factor, following from (!8l- (jT3l . are 

o = l +l \ H ^ [n (3n - 1) B\y\ +C] + C U (51) 
(3n — 1) n 

<P= [ n (3n- l)B\y\+ C]~» , (52) 

where constants C and C\ are not defined by the equations. We will consider only such values 
of parameters n, B and C that the expression [n (3n — 1) B\y\ + C] is positive for any value of 

V- 
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Let us suppose that that "our" brane is that at the point y = (not that at the point 
y = L, as in the previous example). In this case we should take 



Ci 



uj + 1 + n 
(3n — 1) n 



ln(C) 



(53) 



to make four- dimensional coordinates be Galilean on this brane (i.e. a(y = 0) = 0). The 
constant C is defined by the stabilizing potential on the brane at y — 0, whereas the size of the 
extra dimension is defined by the stabilizing potential on the brane at y = L. 

After some algebra we can easily get from (j4"2l) the value of effective four-dimensional Planck 
mass on the brane at y = 



M 2 Pl 



e 2a dy 



C 1 



(54) 



B (2a; + 3 + 3n 2 - 2n) 



2uj + S+3n z -In 

n(3n - 1) BL \ 

c +1 



We will show why such solutions can be interesting from the point of view of hierarchy 
problem by utilizing the choice^ 

n = -. (55) 



Bulk and brane potentials, corresponding to such a choice of n, have the form 

V((f>) = -(3a; + 4)- 



16a; + 5 B3 4 



4 



A 1)2 = T^B (3a; + 4) 02. 
Solutions for the scalar field and the warp factor, following from (151j) . ( 1321) and ( |53l) . are 

4a; + 10 



a 



21 



[In (2L%| +D) -\n(D)\ 



(56) 
(57) 

(58) 
(59) 



21B\y\ + D, 

where D is a constant which will be defined below. 

To stabilize the size of the extra dimension, we add potentials of the form (|38|) on the branes: 



AAl )2 = 71,2 (0 - Ul,2) • 

As in the case discussed in the previous subsection, we get 



\y=0 = Vi, 



\y=L = V 2 . 



(60) 



(61) 



1 The authors are grateful to K. Farakos and P. Pasipoularides for suggesting to examine this case and the 
corresponding background solution, which resulted in this subsection 
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The constant D appears to be denned by 

d = 4( Vl y 3/2 

and we get 

4> = vi 



21(v 1 f/ 2 B\y\ 



+ 1, 

It follows from (IBTj) and fl63|) that the size of extra dimension is defined by the relation 

3 /2 /„, \3/2' 



(62) 



(63) 



215 



(viY /z - fa) 



(flf2 



,3/2 



(64) 



The warp factor has the form 



4q»+10 

= | 21jr : j :i2 /; // ^ J 



(65) 



and it is not difficult to calculate the value of the four- dimensional Planck mass on the brane 
at y = 0. Using (1541) we can easily get 



Mfn = 



<e 2a dy 



B( Vl ) 1/2 (8lu + 27) 



8^+27 



- 1 



(66) 



Let us suppose that all the parameters of the model, made dimensionless with the help of 
a fundamental scale in the TeV range, do not contain a hierarchical difference. For example, 
one can choose B « !TeV~ 7 ^ 2 , v% ~ ITeV 3 and ^ ~ 3.4. It follows from this assumption 

that B (vi) 3 ^ 2 L m 1, which means that L « lTeV~ l . If one chooses ui = 110, then the four- 
dimensional Planck mass on the brane at y = appears to be of the order of Mpi ~ 10 19 GeV. 
Thus, although the largest dimensionless parameter in the model is u = 110, we get a difference 
in 16 orders of magnitude between four- and five-dimensional energy scales. We also see, that 
in this case the hierarchy problem is solved because of the large power in (166]) . contrary to the 
case discussed in the previous subsection, where the hierarchy problem is solved due to the 
exponential factor (as in the Randall- Sundrum model). Of course, if we take larger value of ^, 
the value of uj can be much smaller. 



4 Conclusion and final remarks 

In this paper we considered five-dimensional Brans-Dicke theory as a basis for building different 
solutions corresponding to brane world models. It appeared that for particular choices of the 
potentials and certain values of parameters the theory reproduces some known background 
solutions. We also presented new solutions for stabilized brane worlds, one of which has a 
relationship with a known solution in another frame. We hope that appropriate choice of the 
function F(<f>) can lead to other interesting solutions, which are not evident at the first glance. 
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A reasonable question arises - what happens to the mass of the radion and its coupling 
to matter on the brane, especially in the stabilized cases discussed in subsections 3.3, 3.4? It 
is clear that the radion mass should be expressed through the model parameters 7^2, A (or 
B, respectively), uj and v\,2- Calculations made in [18] for the stabilized brane world model 
proposed in [10] suggest that with an appropriate choice of these parameters the radion mass 
can be made to be in the TeV range, which can be interesting from the experimental point of 
view and does not contradict the known data. Nevertheless, an answer to the question posed 
above can be obtained only after a thorough examination of linearized gravity in the models. 
This issue calls for further investigation. 

We would like to note that all the models found in the present paper have a finite size of 
extra dimension and are of interest for solving the hierarchy problem. At the same time, it 
would be interesting to look for "fat brane" solutions in five-dimensional Brans-Dicke theory, 
analogous to that found in [21], and to examine their properties. This problem also deserves 
to be studied. 

Finally, it is necessary to mention that there are other types of non-minimal coupling of 
the scalar field to gravity. For example, in recent papers [23 [2HJ [2Z] some interesting brane 
world solutions (both analytic and numerical) in the theory with Ricci-coupled bulk scalar field 
were found. In this connection it would be interesting to compare the physical consequences of 
the models with different types of non-minimal coupling and identical solutions for the warp 
factor - for example, of the model discussed in subsection 3.3 of this paper and of the models 
discussed in [HI [23 [21] , in the case of the standard coupling of gravity to matter on the branes 
in all the models. 
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